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Abstract 

Properties of strongly correlated two-dimensional (2D) electron systems in 
solids are studied on the assumption that these systems undergo a phase tran- 
sition, called fermion condensation, whose characteristic feature is flattening 
of the electron spectrum e(p). Unlike the previous models in the present 
study, the decay of single-particle states is properly taken into account. Re- 
markably, the value of the topological charge {N = 1/2) remains unchanged, 
supporting the view that systems with a fermion condensate form a separate 
class of Fermi liquids. Results of calculations are found to be in qualitative 
agreement with ARPES data. 



1 



Typeset using REVT^ 



The single-particle dynamics of Fermi systems at near zero temperatures T is known to 
depend crucially on the index z/, that characterizes the imaginary part of the mass operator, 
\m'L{e 0) ~ e'^ , the energy e being measured from the chemical potential /i. In ordinary 
homogeneous Fermi liquids such as nuclear matter and liquid ^He, where the exclusion 
principle "leads the dance", the index z/ equals 2, and the Fermi liquid can be treated as a 
gas of interacting "immortal" quasiparticles, the cornerstone of standard Fermi liquid theory 
(SFLT) III]. After many successful years, SFLT is currently encountering serious difficulties 
in treating normal states of 2D high- Tc superconductors. The analysis of ARPES data 
shows that even around the diagonals of the Brillouin zone, the index v is unity @J3, while 
in the immediate vicinity of the van Hove points (vHP), the sharp ARPES peaks disappear 
altogether [|^. 

We propose that solution of this challenging problem is associated with fermion conden- 
sation [[^-[T7[|, a novel phase transition that generates a group of dispersionless states, called 
the fermion condensate (FC), whose energies coincide with /x. States with the FC have been 
uncovered as unconventional solutions of equations of Landau theory at T = 0. To gain 
insight into the problem, consider the Landau formula for the variation of the ground state 
energy Eq 

SEo = J2 e(P; ^f) 5n(p) + ^ H r(pi. Pa; np) 5n{pi) 5n{p2) , (1) 

P Pl,P2 

applicable if the change of the quasiparticle momentum distribution ripip) = 6{fi — e(p)) is 
rather small. The single-particle (SP) spectrum e(p) and the scattering amplitude r(pi, pa) 
are to be evaluated for the Landau state. The formula (1) is appropriate for homogeneous 
systems, where p is the momentum, as well as for finite ones, where p stands for a set of 
relevant SP quantum numbers. The crucial point is that due to the discreteness of the SP 
spectra, low- lying states of finite systems, adjacent to the ground state, do not practically 
damp. This situation persists until the excitation energy e attains values ~ e%/^/A, (e^p is 
the Fermi energy, A is the particle number), above which the SP widths 7/(e), calculated 
with incorporation of broadening of energy 5-functions, exceed the distance between the SP 
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levels. Then the widths 7/(5) and 7(e), inherent in infinite systems, rapidly approaches 
each other. We infer that in search for a new ground state of a large system, the momentum 
representation can be employed, while damping effects can be ignored. At the same time, in 
dealing with the spectral functions at realistic energies, damping effects should be properly 
taken into account. 

As seen from Eq.(|l]), the Landau state holds until the necessary stability condition (NSC), 
requiring the positivity of the l.h.s. of Eq. (|l]) at any admissible variation of np-, is violated, 
and thereby the minimum of Eo{n) moves away from np- A prerequisite for this to occur is 
the presence of substantial momentum-dependent components in F (see below). For a long 
time, this problem was studied only under the assumption that the new Fermi surface be- 



comes multi-connected, while the occupation numbers n(p) remain 1 or [|T8|. However, the 
expansion (|l|) can also be applied to evaluate variations of Eq when continuous quasiparticle 
distributions n{p) come into play, and beyond the point where the NCS is violated, it often 
turns out that such unusual for SFLT, but accustomed e.g. for the Bogoliubov theory of 
superconductivity, distributions yield a lower energy than any discrete-valued solution with 
multi-connected Fermi surface. This compels us to treat the quasiparticle occupation num- 
bers n(p) in Eq. (|I|) as variational parameters and look for a new location of the minimum 
of E{n) from the variational condition 

SEojn) 

jr I X = , P e C , (2) 
(5n(p) 

with 5Eq given by Eq. (|l]). The region C is the place of the residence of the fermion 
condensate, since the l.h.s. of Eq. (0) is just the quasiparticle energy e(p; n), and Eq. (0) is 
equivalent to e(p G C; ra) = yU. Outside this domain, whose boundaries are determined from 
Eq. (^ itself, the old solution npijp) survives. 

Similar results can be obtained by summation of diagrams of the mass operator in Dyson 



equation The transition occurs, if in this equation, there exists a bifurcation point. 



Beyond this point, the ground state distribution n(p) is evaluated from relation 

e°+S(p,£ = 0;n)=/i, peC, (3) 



which is equivalent to the equation e(p G C;n) = fi. If one neglects for a while the e- 

dependence of E(p, e), then the FC Green function G(p G C,e) = — 6° —S(p G C, e)]^^ 

becomes 1/e, that leads to a change of the Fermi liquid topological charge N 0. Recall that 

for ordinary homogeneous Fermi liquids, with G{p,e) ~ [e — e{p) + /i]"^, and for Luttinger 

liquids as well, the value of is unity, as given by the integral 

r dl 

N=i— G(p, e=in) diG-\^, e=in) , (4) 
J Zm 

L 

taken along the contour L, embracing the Fermi line in the 3D space {pxiPy-,^)- However, 

upon inserting Gc{p, s) = i/£ into the integral (H), one finds = 1/2, implying that systems 

with the FC form a separate class of normal Fermi liquids As we shall see, this result 

holds if the energy dependence of the mass operator is incorporated. 

The substantial momentum dependence of F, needed for solution of Eqs. (H) and @ 

to exist, emerges, for example, in the vicinity of an impending antiferromagnetic phase 

transition. In this case, the static spin susceptibility x(q, = 0) ~ [P^ + /t^(q — Q)^]~^ 

||T9| diverges at q — > Q = (tt, tt), since (3 vanishes at the transition point. The direct term 

Frf in F, being proportional to — x(q, = 0), diverges as well, and so does the exchange 

term Fe(p, pi, q) ~ [P^ + fi;^(pi — P2 + q — Q)^]^^ because of the antisymmetry relations 

imposed on F [^. It is the exchange term, taken at q = 0, that enters Eq. (|T]). Parameters 

/3 and k are assumed to depend on the density p rather than on ^^^(p) itself, so that 

6Te/Snp{p) = 0. If the amplitude F is approximated by Fg, then upon inserting Eq. (|l]) 

into Eq. (0) the latter is recast to ei(p; np) +J2^e{p, Pi; np) n{pi) = fi. From the definition 

pi 

ei(p) = e(p; ^f)— re(p, Pi)?t-f(Pi)5 it is clear that the variational derivative 6ei{p) / 6np{p) 
pi 

vanishes. Therefore ei(p) is reduced to the LDA spectrum ep, and Eq. (|^) takes the form 

4 + Z!re(p,Pi)r^(pi) =/i , peC (5) 
PI 

often employed in theory of fermion condensation. It is worth noting that Eq. can be de- 
rived from the effective energy functional Ees{n) = J2^p n{p) + 1/2 J2 ^e{p, Pi) ^(p) ^(Pi)j 

p ppi 

which is appropriate only for the rearrangement problem and has no connection with the 
Hartree-Fock part of the ground-state energy functional. 
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Results of numerical calculations of Eq. (^), which are insensitive to the choice of the 
parameter k, show that a FC arises when the filling approaches 1/2. Once it appears, the FC 
resides close to the vHP. As seen from Fig. 1, the fraction rj of the Brillouin zone occupied by 
the FC remains rather small, attaining a maximum ~ 0.1 when the filling slightly exceeds 
1/2, and the hole pocket is centered around (vr, tt) (for more details, see IT^). 



The degeneracy of the SP spectrum at T = 0, a salient feature of the solution given by 
Eqs. (H), (H), is lifted by pairing interactions which are not included into Eq. (P 
doing so, the BCS occupation numbers f^(p) coincide with n(p) evaluated from Eq. @ or 
Eq. (0) provided the BCS coupling constant is rather small. Evidently, in obtaining both 
these distributions damping effects can be ignored. For this reason, superfiuid systems with 
and without the FC look more alike than normal ones, since in normal states of conventional 
Fermi liquids the damping makes no difference, whereas in normal states of systems with 
the FC, the damping becomes a real "weathermaker" . Indeed, the relevant contribution to 
ImSij(p, e) is given by ||21[ 



q.pi 



ImS/j(p,£:) ~ ^ // dujdeiF{e,uj,ei,T)\r{p,e,pi,ei,ci,uj;n)['' 



X ImGn^p - q,e - uj)lmGR{-pi,-ei)lmGR{q-pi,iJ - El) , (6) 

where the factor F{e,uj,ei,T) = cosh(^)[cosh(|i;) cosh(^^) cosh(^i^2^)]"^, |rp is the ab- 
solute square of the scattering amplitude, properly averaged over spin variables, and Or is 
the retarded Green function with ImG/j(p, e) = — 7(p, e) [{e — cr(p, e) — s(p))^ + 7^(p, £)]~^, 
where a{p, e) = Re S_r(p, s) — Re Si?(p, e = 0) and s(p) = Cp + Re Sij(p, e = 0) — fi. 

We restrict ourselves to temperatures T markedly lower than the characteristic temper- 
ature Tf of destroying the FC that allows one to ignore the T-dependence of rj. To simplify 
the problem, we replace the function 7(p, e) by a set of functions of the single variable e, e.g. 
in the FC region, 7(p G C, e) is reduced to 7c(^)- The complementary region of momentum 
space, in which the dispersion of the spectrum e(p) has a nonzero value, is composed of two 
subregions. The first, adjacent to the FC domain and henceforth denoted by T, is a transi- 
tion region, in which the same decay processes, as in the FC region, are still kinematically 



allowed. The second subregion, denoted by M, is located around diagonals of the Brillouin 
zone. Here some of these processes are either kinematically forbidden or at least strongly 
suppressed. Correspondingly, 7(p G T,£:) — > 7t(£^),7(p G M, £:) — > •yui^)- To close the set 
of equations of the problem, the amplitude F should somehow be specified. Bearing in mind 
that 7] is small, we shall initially neglect the FC contributions to F, replacing it by F(ni?). 

We start with the case e ^ T, and set T = in the integral ®, thus dropping all T- 
dependent contributions. First we evaluate 7c (^^ ~* 0)- this case, (i) contributions from 
processes involving only the FC states prevail (see below), (ii) the quantity |F(ni;')p ~ 
can be factored out of the integral @, and (iii) the quantity s(p G C,T), which vanishes 
over the FC region at T = 0, can be verified to remain small compared to leading terms in 
Eq. (P), and thus can be neglected. As a result, the energy and momentum integrations in 
Eq. d^) separate. Taking for certainty e > and omitting numerical factors, we are left with 

7c(£^0) ~/3-V / / Ac{e-uj)Ac{-e,)Ac{uj-e,)du;dei , (7) 





where Ac{e) = ImGij(p G C,e). To proceed, we insert 7c(^ ^ 0) ~ e^*^ into the Kramers- 
Kronig relation to obtain ac{£ ^ 0) ~ e'^'^. We then substitute 7c and ac into Aq and find 
Ac{b ^ 0) ~ e'"^. Finally, upon inserting this result into Eq. (0), we arrive at uq = 1/2 
|T5| . More precisely, one obtains 7c ^ 0) ~ /^^^(r^e^e)^/^ and 

G«(p G C, £ ^ 0) ~ e-'^'/^Ms ^ 0)]-^ ~ e-'^'/^Pive^er'/' . (8) 

This result can be shown to hold even if the momentum dependence of the quantities 7c (p, £) 
and cTc(p,e) is properly taken into account. We see that in the FC region, the conventional 
structure of the Green function is destroyed, the familiar pole being replaced by a branch 
point at e = 0. What happens to the topological charge A^? Upon inserting the Green 
function into Eq. and performing simple integration, we are again led to Volovik's 
previous result N = 1/2 in spite of the dramatic alteration of the Green function itself 
that occurs in the FC region. 

In the transition region T, the decay into the FC states is not kinematically forbidden. 
Accordingly, ^^{e 0) ~ {r]e'^pey^'^ , while the function s(p G T) already differs from 
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zero. Requiring it to vanish at the boundaries of the FC region along with its first derivative, 
one finds that in the region T, the conventional structure of the Green function is recovered, 
but in the vicinity of the FC domain, single-particle excitations appear to be ill-defined, 
since the pole of G{p, e) is located close to the imaginary energy axis. 

In the region M, dominant contributions to Im come from a process associated 

with the generation of three states: two from the FC region and one from the M region. In 
this case, the formula for finding 7m (^^ — ^ 0) reads 

£ U) 

7M(£^0)~r'E / I Ac{-ei) Ac{uj-ei)[l - e{Y*)]Pc{v-Vi) Au{v,e-uj)dudei , (9) 



p,pi 







where -Pc(ct) = Z]p^(p)^(p ~ q) ^-iid 0(p) = 1 if p G C and otherwise vanishes. It is seen 
that in this case, the momentum and energy integrations do not separate. However, one can 
take advantage of the fact that the spectrum ^m(p) = ^m{v) — is proportional to {p — Pf) 
and introduce ^m(p) as a new variable. Then after simple integration, we are led to the 
result i^M = 1 postulated in the model of a marginal Fermi liquid (MFL) p2[. Evaluation 
of the ?7-dependence of relevant quantities in the M region yields 7m(£^ ^ 0) ~ (3~'^r]^^'^e and 

These results can be applied to the case e ~ T, where according to Eq. (^, the leading 
term in the FC Green function has the form Gr(p G C,e) ~ e~'"'/^[7c(£^, Upon 
inserting this expression into Eq. (|^ one finds that the damping 7c(x,T), where x = e/T, 



can be displayed as 7c (a^, T) = '-fcyTsp D{x) where the constant 7c specifies the compound, 

while the dimensionless quantity D{x) obeys the integral equation 

, - _ v,^ /"/" dydz , , 

(a;j = cos -JJ ^^^^ K cosh ^ cosh ^D{-y)D{x-z)D{z-y) ' ^ ^ 

With this result, the damping 7M(a;, T) is calculated straifghtforwardly: 

, . X [ f dydz , . 

7M(a:,T) = 7MTcosh-yy f cosh ^ cosh ^ ' ^^'^ 

the constant 7m being a characteristics of the given compound. The function 7M(a;, T)/T 
starts out of the origin as a parabolic function 7m(x, T)/T ~ 1 + 0.1 x^. The asymptotic 
regime 7M(a;, T)/T ~ x, stemming from Eq. (|), begins at x ~ 2.5. 



In normal states with the nonzero value of the pseudogap, relations (|TOD, (|TT|) hold, as 
long as it remains less than the damping 7c (T). These relations persist even in superfiuid 
states provided the gap value meets the same restriction. On the other hand, they are 
violated if energy attains values, at which contributions to 7(e) that were omitted from 
Eqs. (0) and (|^) become comparable to the terms that were retained. A leading correction 
to the integral comes from final states, that involve one hole (particle) belonging 
to the region T. Eq. (P) can be employed to estimate this contribution, with the single 
replacement s(p G M) s{p G T). We find ^ 0) ~ P^'^s, which is independent 

of the 7] value. Estimation of other corrections to 7(e) is carried out along the same lines, 
justifying the identification of (|^) and (H) as paramount contributions to lmT,^{p,e — > 0) 
until e exceeds the characteristic PC energy epc — V^'py evaluated by comparison of S'jci^) 
and 7c(£) ~ f^-'ivep^y/^ 

At energies e > epc, the corrections exhibit themselves in full force, so Eq. should be 
solved numerically in conjunction with the Kramers- Kronig relation, employed to connect 
7(e) and cr(e). This is done with the aid of an iteration procedure, which converges rapidly. 
In Fig. 2 we display results from these calculations carried out with the same amplitude F 
that was employed in the calculations of the FC characteristics shown in Fig. 1. Two different 
rj values, specifying the fraction of the Brillouin zone occupied by the FC, were considered: 
(a) 1] = 0.1, close to the maximum rj value in the model of fermion condensation driven 



by antiferromagnetic fluctuations |jT^, and (b) rj = 0.01. In spite of the simplicity of the 



interaction adopted, salient features of ARPES data ||3|,|6|JT^ are reproduced, including the 
universal behavior of the ratio Im T,^{p G M, e, T; 77) /T as a function of x = e/T, uncovered 
in 0. Moreover, our theory predicts the same behavior of lmJ]ji{p G M,e,T;ri)/T for 
different compounds provided results are properly normalized. We choose this normalization 
to ensure the same slope k{ri) = |9ImS/j(p G M,e,T]ri)/d6\ at a; ^ 1 where the damping 
changes linearly with e. The latter is demonstrated in Pig. 3, where two functions |Im Sij(p G 
M, e, T; ri)\/T, evaluated at 77 = 0.1 and r] = 0.01 and normalized to the same slope k = 0.75, 
are compared to data measured in |^ for the optimally doped cuprate Bi2Sr2CaCu208+5. At 
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the same time as follows from Eq. (|TOD, in the FC region the above universal scaling of the 
ratio lmTjji{e,T)/T is destroyed, and instead ImSij(p G C,e,T) displays -\/r-dependence 
at a; < 1. 

The above scenario in which the fermion condensation precedes the antiferromagnetic 
phase transition does apply in the three-dimensional case, although the range of the FC 
region shrinks markedly. Along the same lines, one can analyze the situation with fermion 
condensation in the vicinity of other second order phase transitions, such as charge-density- 
wave instability. So far the feedback of the FC on the scattering amplitude F has been 
ignored. However, the simplest FC diagram, i.e. a loop, evaluated with the FC Green 
function (H), diverges logarithmically. As a result, we are led to a familiar problem of the 
parquet- diagram summation, solution of which will be reported in a separate paper. 

Summing up the results of our analysis, we infer that infinite electron systems with a 
fermion condensate, independently of the dimensionality, do not admit Landau quasiparti- 
cles, since the renormalization factor z = (1 — dT,/de)p^ that determines the quasiparticle 
weight in the single-particle state vanishes in all the regions of the Brillouin zone. In the FC 
domain, the value of the topological charge N is found to be = 1/2 indicating that sys- 
tems with a fermion condensate form a separate class of Fermi liquids. The model of fermion 
condensation presented here allows one to explain basic features of the spectral functions of 
normal states of high-Tc superconductors, including the MFL behavior of the damping of SP 
states around the diagonals of the Brillouin zone (the M region) and the suppression of the 
peaks in APRES data in the immediate vicinity of the van Hove points. And the universal 
behavior of the ratio ImE/j(p G M,x,T)/T as a function of x = e/T established in for 
the M region is also reproduced in this model. Moreover, when properly normalized, our 
results for ImS^(p G M,x,T)/T calculated for different compounds collapse on the same 
curve. However, as we have seen, this universal behavior is destroyed in the C region close 
to the van Hove points. 
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FIGURES 




FIG. 1. The Fermi surface in the model of fermion condensation, driven by antiferromagnetic 
fluctuations with the scattering amplitude r(p,pi) = (Af(0))~^[/3^ + K^(p — Pi — Q)^]~^, where 
N{0) is the density of states and /3 = 0.2. The regions occupied by FC are shaded. 
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FIG. 2. The damping of the single-particle states at T = in the vicinity of the vHP 
(solid lines) and around the diagonals of the Brillouin zone (thin lines), calculated with the same 
scattering amplitude F for 77 = 0.1 (upper panel) and 77 = 0.01 (lower panel), and measured in s'p. 
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klk{r]) ■ |ImEjj(peM,e)|/T 




s/T 



FIG. 3. The ratio |Im S/j(e)|/T around the diagonal of the Brillouin zone as a function of 
e/r, calculated for r] = 0.1 (solid line) r] = 0.01 (thin line) and normalized to the same slope 
k = 0.75, i.e. multiplied by the factor k/k{rf). The experimental data for the optimally doped 
cuprate Bi2Sr2CaCu208+5 @ are shown by open and solid circles (T = 300 K), triangles and 
squares (T = 90 K), and diamonds and stars (T = 48 K). 
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